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Abstract 

The study of Deeply Virtual Compton Scattering has shown that electromagnetic gauge invari- 
ance requires, to leading order, not only twist two but additional twist three contributions. We 
apply this analysis and, using the EUis-Furmanski-Petronzio factorization scheme, compute the 
single (electron) spin asymmetry arising in the collision of longitudinally polarized electrons with 
hadrons up to twist 3 accuracy. In order to simplify the kinematics we restrict the actual calculation 
to pions in the chiral limit. The process is described in terms of the generalized parton distribution 
functions which we obtain within a bag model framework. 

1 Introduction. 

Hard reactions provide important information for unveiling the structure of hadrons. The large vir- 
tuality, Q^, involved in these processes allows the factorization of the hard (perturbative) and soft 
(nonperturbative) contributions in their amplitudes. Therefore these reactions are receiving great at- 
tention by the hadronic physics community. Among the hard processes one, which merits to be singled 
out, is the Deeply Virtual Compton Scattering (DVCS) because it can be expressed, in the asymptotic 
regime, in terms of the so called Generalized Parton Distributions (GPDs) [||, |3|. The GPDs de- 
scribe non-forward matrix elements of light-cone operators and therefore measure the response of the 
internal structure of the hadrons to the probes. Moreover DVSC is instrumental in the experimental 
interpretation of the angular momentum sum rule Q . 

It has been shown that the implementation of gauge invariance in the analysis of the DVCS 
amplitude in the asymptotic regime, i.e. large virtuality of the incoming photon, requires the inclusion 
of twist-3 contributions ||^]-|jl3[. Let us explain the reason in a brief manner. In leading twist, in the 
Bjorken limit, the Lorentz structure of the hard subgraph of the DVCS amplitude's leading diagram 
has, at large Q^, the form of a transverse projector. The virtual photon momentum in the form of the 



Sudakov decomposition contains a transverse component too. Thus their contraction does not vanish 
and electromagnetic gauge invariance is violated. To restore it, next-to-leading order terms in the 
asymptotic expansion to the DVCS amplitude, which are proportional to the transverse component of 
the momentum transfer, have to be included. These terms, which are twist-3, give rise to the dominant 
contribution in some observables. One of them is the single spin asymmetry (SSA), which arises in 
the collision of longitudinally polarized electrons with hadrons, and which we will analyze in here 

In this first paper we deal with a spin zero massless target, the pion. The same procedure could be 
applied to the scattering off unpolarized nucleons. However this calculation would require keeping the 



mass terms of the nucleon, a complication which we want to avoid at present [16|. Moreover, in the case 
of polarized nucleons, one could study besides the SSA other asymmetries, however the complexity of 
the analysis, with the existence of many different GPDs, is postponed for a future publication [|16|. In 
order to calculate the twist-2 and twist-30 GPDs contributing to the SSA, a crucial ingredient of the 
calculation, we use an MIT bag model scheme with boosted wave functions. 

The plan of our paper is as follows. In section 2 we describe the kinematics and introduce the 
appropriate notations. The DVCS amplitude for the pion, including up to twist-3 contributions, is 
presented in section 3. We follow the description of Ref.0], where, using a generalization of the Ellis- 



Furmanski-Petronzio (EFP) factorization scheme |17], the complete gauge invariant DVCS amplitude 



has been obtained. In section 4 we outline the basic ingredients of our approach and calculate all 
the parametrizing functions (GPDs). Then in section 5 we give the numerical estimates of the SSA 
parameter for the case under study and discuss our results. 

2 eh ^ ^eh process: kinematics and notations. 

Our starting point is electron-hadron scattering into real photon, electron and hadron, 

e(A:)+vr(p) ^7(^') + e(A:')+vr(p'), (1) 
assuming that the electron is longitudinally polarized. We will consider the electron to be massless, 

= k'^ = 0, (2) 
and, that the hadron is a pion, which we take also to be massless (chiral limit). 

In QED, to lowest order, the reaction (|^) takes place via the Bethe-Heitler process, where the final 
real photon is emitted by one of the electrons, and the virtual Compton process 

7*(g)+7r(p) ^7((7')+7r(p'), (4) 

which becomes the DVCS process if the square of virtual photon momentum q = k — k' is very large, 
i.e. 

^Recently, a detailed report with the description of all possible asymmetries has appeared in the literature (see 
Ref.0]). 

^We stress that throughout this paper we mean twist-3, which include both the kinematical and dynamical twist-3 
contributions. 



2 



For the Bethe-Heitler process the ampHtude is given by 

= -e'^(^0{^*(|f^7/. + 7,<(^^;^e*}«(A;)^r^(p,p'), (6) 

where 

A = p'-p, (7) 

and the electromagnetic vertex of the pions takes the standard form 

T^{p,p') = {p + p')f,F+{A^). (8) 

i<+(A^) is the electromagnetic form factor of the pion. The Virtual Compton process amplitude is 
given by 

Tivc) = ^(fc')7^^(A;)^T/^/^^e;, (9) 

where T^^'^'^ corresponds to the DVCS subprocess. 

To describe the reaction (^) it is useful to introduce the following dimensionless fractions 

p-q 

X = , y = — r, 

2p ■ q p ■ k 

p-p' p-q' , . 

z = or . (lUj 

p-q p-q 

These fractions can be related with the Mandelstam variables of reaction (|l|) and of the DVCS sub- 
process. Indeed, if we introduce the following variables 

s = {q + P?, i={p'-pf, (11) 

for the DVCS subprocess, and 

S = {k + pf ^2k-p, (12) 

for the reaction (||), the following relations hold 

= xyS, s = {1 — x)yS, i = ySz. (13) 

Since we neglect the pion mass, the most suitable system of reference is the center of mass system, 
where 

A; = £;(l,sin/?,0,cos/3), q= {Qo,0,0, -Ep), 

p = Ep{l, 0,0,1), p' = Epi{l,sm6cos(f),sm9sm(j),cos6), (14) 

and 

Ep='^^^, Ep,=E,, = ^V~s, (15) 

Moreover, we define (p to be the angle between the plane formed by the three-dimensional vectors 
q' and k (leptonic plane) and the plane formed by the three-dimensional vectors q' and p' (hadronic 
plane). It can then be easily seen that the other angles satisfy the following relations 

cos6' = l-2z, cos/? = 1 - ^. (17) 

1-xy 
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3 DVCS amplitude off pions up to twist-3. 



In this section we focus on the DVCS amphtude off pions. We start from the expression for the virtual 
Compton scattering amplitude which can be written as usually in the form 

Tf,u=iJ d'^ze'i'' {p'\TJ^{0)Mz)\p), (18) 

where is the electromagnetic quark current: 

Jf,{x) = i^{x)Qj^i^{x). (19) 

In Eq.(|l^) Q is the charge quark matrix, which is equal to 

^(l + Srs) for5C/i.(2) (20) 

and to 

^(Aa + ^As) for SUf{^). (21) 

Here Tj and Aj are the conventional Pauli and Gell-Mann matrices for two and three flavors respectively. 

As discussed before a gauge invariant DVCS amplitude cannot be written down unless the twist 
three contributions to the amplitude are taken into account. Here we would like to recall, briefly, the 
results obtained in Ref.p and reproduced by several groups 0, ^. 

The main point of these analyses is that, the violation of the photon gauge invariance, is propor- 
tional to the non-zero transverse component of the virtual photon momentum |^. In other words, 
the convolution of the leading order DVCS amplitude with the virtual photon momentum is propor- 
tional to the first degree of transversity that corresponds to twist-3 |^. Hence to obtain the complete 
gauge invariant DVCS amplitude we have to consider all terms which are linear combinations of the 
transversity. This situation is not surprising and similarity with the transverse polarization in deep 
inelastic scattering off nucleons can be recalled [19, [20|]. 



Therefore, in order to preserve the electromagnetic gauge invariance up to leading order within the 
generalized EFP factorization scheme, one must add to diagram (a) of Fig.||, consisting of a hard part 
with two quark legs the diagram consisting of a hard part with two quark legs and one transverse 
gluon (see diagram (b) of the same Figure). This latter diagram is entirely twist 3, while diagram (a) 
contains, besides the standard twist-2 term produced by the good components of the quark fields and 
collinear parton momenta, a twist-3 term, which can be related to the quark gluon contribution of (b) 
by means of the equations of motion. 

After performing the T— product for the two electromagnetic currents and going from the four- 
dimensional integration over z to the one-dimensional integration over the x— fraction the amplitudes 
of diagrams (a) and (b) may be written as |^, 18|: 



J dxidx2trl^Eij,p,y{xiP,X2P)ujpp'Tp>{xi,X2)^ , (22) 



^In the present case, the transversity is the transverse component of momentum transfer. 

^In this paper we deal with the Born diagrams only. The EFP factorizing scheme for the general case can be found 
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where ujppi = 5ppi — Up/Pp, and 



Epu{xP) = jpS{xP -j + q)lu + ''crossed' 

A A 

Epp^{xiP,X2P) = -fpS{xiP - y + q)lpS{x2P - y + 9)71^ + "crossed" 

Nf 



iV f 

T{x) = - [ dAe*(--+«)^ J2 {p'\M>^n)MO)\p), 

•' a=l 

TP'{xi,X2) = -J dXidX2 e^(^'i+€)^i+^(^2-xi)A2 

J2 el {p'\M>^in) DP' (A2n)4(0)|p). (23) 



a=l 

Here Dp is the QCD covariant derivative in the fundamental representation ^. 

The use of the QCD equations of motion lead to the following expectation values 



{p'\D{z)i^{z)i>mp) = 0, (pXz)V;(0)D(0)|p) = 0. (24) 

Keeping only the non-zero (axial and vector) projections of the quark and quark-gluon correlators, 
we are able to express the tree-body (quark-gluon) parametrizing functions in terms of the two-body 
(quark) parametrizing functions. As a result of this trick, the amplitude corresponding to diagram 
(b), expressed by means of two-body parametrizing functions, grouped together with the amplitude 
corresponding to diagram (a), expressed also by two-body parametrizing functions, leads to a gauge 
invariant DVCS amplitude, which reads 



rpDVCS I f ^ I , I W 

~ 2P-qJ [x-^ + ie^ x + ^-ieP"" 



(25) 



where 



Tp, = Hi{x) \^-2CPpP, - PpQ, - PuQp+ 

5^.(P-Q)-^pX + ^^-Am) - 

H,{x) (^^P.Aj; + S^PpA'^ + AlQ, + AlQp^ - 

^Ha{x) (^SeP^A^ - ^P.Aj; - AlQ, + AlQ^j . 
Here the parametrizing functions (GPDs) are, 

Nf 

H,{x) = elH^{x), i = {1, 3, A}, Nf = number of flavors, (26) 

a=l 

which were defined in Ref.|^], and whose calculation we show in Section 5. The DVCS amplitude is 
complex, due to the factor in front of Tpi, in Eq.(^), leading therefore to a non vanishing SSA, which 
could be seen in the collision of longitudinally polarized electron beams with pions. 
^We follow the notation of ref.Hl 
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4 Single spin asymmetry parameter. 

We next calculate the SSA parameter |p as a function of the angle (f) between the leptonic and hadronic 
planes. The SSA parameter is given by 

_ d(j{^) - da{^) 

where d(T(— ) denotes the differential cross section with different helicities for electrons. 

The difference of cross sections in the numerator of (^) is defined by the imaginary part of the 
convolution of the leptonic tensor with the hadronic tensor and consists of two terms 

da{-^) - da{^) = Ada^^\^; ^) + Ada^^\^; ^). (28) 

The first term in (^), Ada^^\^; emanates from the interference between the Bethe-Heitler and 
the virtual Compton processes. Its contribution is equal to 

Ada^'\^; = idPSf2-^MX^^hmHg-i^^\ (29) 

where {dPS)^ is the three-particles phase space, and the leptonic tensor L^j^^^^ is given by the 
following trace 

^S"^ = tr (7m^{7. (^^7« + 7« (|f^7.}^') ■ (30) 

Finally, the hadronic tensor H^J^a^^ is given by 

HgT^ = T^/'^'ip + p')aF+{t). (31) 

The second term in (^), Ada^^\^; <— ), is related to the square of the virtual Compton amplitude 
and is defined by the expression 



Ada(^)(^; ^) = idPSf2^L(fJlmHjfJ, (32) 



(s) 

where the leptonic tensor Ll^u is 



4^.^=tr 7^75^7.^:' , (33) 



(S) 

and the hadronic tensor H^,y is given by 



Hi'J = r,^/^^(T.^/^^)+. (34) 



Calculating the traces in ( |30D and (^3|) , and the imaginary parts of (|3l| ) and (^^ , which arise in the 
DVCS amplitude as discussed before, we obtain for the first term in the difference of cross sections |^] 



J dx(^6{x + ^)-6{x-C)j • (^Hi{x){{k + k')- P) + 



2H,{x){k' ■ A^) + -^^^Ha{x) ( {k ■ A){k' ■P)-{k'.A){k.P)]], (35) 
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and for the second term 



V 



Ada(-'^\^-^) = idPSf^ekk'PA 

6{x + 0- S{x - 
Six' + e) - Six' - e) 



2e 



+ 



iP-Q) 

V 



dxdx' 



x' + C 



V V 

+ 



X — ^ + C 
Hxix)B^ix') - Hxix)H?X^)^ 



' R,ix')^ - H,ix)^^ 



(36) 



Let us recall the notation for (^) and (^): is the pion electromagnetic form factor, arising 

from the Bethe-Heitler diagrams, k and /c' denote the momenta of the initial and final electron. 

A careful analysis of Eqs. p5| ) and (|3^) shows that the twist-3 parametrizing functions (GPDs) 
and Ha appear in ( [35| ) as corrections to the twist-2 parametrizing function i/i, while in Eq.(|36|), 

the twist-3 parametrizing functions, appear in the leading terms, and therefore, it is important to 

emphasize that they give the main contribution. 

Next, we turn to the consideration of the DVCS process off pions with unpolarized leptons. The 

cross section for this case is 



(37) 



dounv = idPSy [\TBHr + \Tvcr + iTBHT^c + Tvcnu) 
The pure Bethe-Heitler contribution to the cross section reads (see pl) 



(38) 



where the leptonic tensor L^^^^ is defined by 

k' + A 



-'flU 

k- A 



k' + A k-A ^^yj 

ik - Af^^ + ^^ik' + Af' rVik' + Af^" + ^"ik-Af' 



and the hadronic tensor H^u^^ is defined as 

i^jr) = r^r, = ip + p'),ip + p').Flit). 

For the virtual Compton process we have 



Q 



,4. fiu -^-^ 111/ 



(39) 



(40) 



(41) 



where the leptonic tensor is given by 



4^^^ = tr 



liikluk' 



(42) 
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and the hadronic tensor is given by 



Here, as usual, T^^*^"^ denotes the DVCS ampUtude. 

The contribution arising from the interference between the Bethe-Heitler and virtual Compton 
processes is given by 

TycT+^ + TbhT+c = (44) 
where the expressions for Ljji!^a'^^ and H^ij^a^^ have been presented above. 



Furthermore, using Eqs.(pq) and ( [37| ) and inserting for the parametrizing functions (GPDs) Hi, 
and their values calculated within the MIT bag model, as will be presented in the next section, 
we obtain, as a function of (f>, the values for the SSA parameter shown in Figs.§, ^, and|^. Note that we 
have chosen the kinematics for the reaction (|^) which can be achieved at the HERMES experiments 



(see also |21, p2| 



5 GPDs within an MIT bag model scheme. 

The Generalized Parton Distributions (GPDs) are defined as the parametrizing functions of light-cone 
matrix elements of bilocal field operators f\ 

J ^e^'''{<p'm-^)m^)w{p)), (45) 

where x is the momentum fraction of the parton, a light-cone vector to be specified later and T 
represents Dirac matrix structures. The functions that we need arise from the expectation values of 7^ 
and 7^75. With the notation of Ref.Q the parametrizing functions appearing in the SSA parameter 
are defined by the following matrix elements 

E^lJ ^e^''{<P')\M-^hM^)\^(j^)) = Hi{x)P, + H,{x)Aj^, (46) 

a=l 
Nf 

E^lJ ^e^''(^(p')lV^a(-^)7M75^a(^)K(p)) = zHA{x)e,.p.P''n^Al, (47) 

a=l 

where ^ 

P = , A = p' - p , = A - (A • n). (48) 

All the parametrizing functions are functions not only of x but also of the momentum transfer 
t = A^ and therefore connect the parton distributions and the form factors ^, |3|. Moreover each 
parametrizing functions possesses the following properties (see for instance Q): if x belongs to the 
interval [^,1] then the //j— functions for fixed flavor are with the quark distributions; if x belongs 
to the interval [— 1,— then the i/j— functions for flxed flavor are the anti-quark distributions; the 
®In eq. the triplet of pion fields is defined by tt = {tv" ,n'^ ,tt^} (cf. 
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iJj— functions are the difference between the quark and anti-quark distributions when x belong to 
interval ] — ^, i.e. 



^ 

a=l 



Nf 



Y.ei^nux,m{-^ <x< 1) - wf(x,e)e(-i <x< (49) 

Here the Tif and Tif related to the b~ and combinations of creation and annihilation opera- 

tors, respectively. 

We next proceed to calculate the parametrizing functions within the MIT bag model. We choose 
the kinematical variables in the Breit frame which become 



p'^ = (M, ^),p^ = (M, , = (0, Ax, -2eM). (50) 



and the light-cone vector is given by 

1 



n^ = =(1,0,0,-1). (51) 
^From all these equations it is easy to obtain the expressions for the parametrizing functions, 

a=l 
Nf 

Hsix,^,t) = ^e2-L-|^e^^-(7r(p')fc(-^)4±V'a(^)|vr(p)), (53) 



a=l 
Nf 



a=l I -'-I 



^e^^H^{p')\^t{-^)^yM^)\<p)), (54) 



where is the y component of the spin operator. We need to emphasize, that within the most 
naive version MIT bag model, i.e. when only confinement is taken into account and no evolution is 
considered, only the valence quarks degrees of freedom are considered. As a consequence, in the Tr"*" 
case the matrix elements of Eqs.(^2|), ( ^3|) and (|5^ reads 



Hf\xX) = ^wr(x,oe(-c <x<i)- hnf{x,i)@{-i <x<i) = 

^ e^^^ {.Hp'm-^) r. V'(^)k+(P)>, (55) 

where Tj = {yi, ^x,7oSy} for the corresponding parametrizing functions. Similar expressions can be 
easily written down for the other terms of the pion triplet. 

In order to perform the calculation we use the MIT bag model in the boosted scheme p^ , whose 
virtues and defects for this type of physics have been thouroughly discussed . The expressions for 
the parametrizing functions above become in this framework of the form 

2M j ^e^'^ J d^xe^^- W- ^) T ^)) (56) 
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where now the matrix elements are calculated within the bag states normalized to 1 and T symbols 
the required Dirac operators. These field operators give rise to a sum over quark (anti-quark) wave 
functions of the form 



2MZ ' 



2tt J 



An, 



An, 



(57) 



where Z is a normalization factor coming from the spectator particle and the boosted wave functions 
are given by 



(2^ 



exp {—Hot — k ■ r}ip{k) 



(58) 



All the required definitions and notation from now on are to be found in Refs.[p3] ^-^d |24]. 

After a tedious but straightforward calculation we obtain the parametrizing functions for each 
quark flavor 



A-kNqR^ZM 



cosha;(l — (cosh a; — 1) 



dk±dip 
"(2^ 



k' -k 



k'k 

+(coshu; - 1) (^j^{vzk- v - kzv'^)(fo{k')(pi{k) 

+77^(^2^' ■ ^ - k'zV^)(po{k)(pi{k') 

f k ^ k' ^ 

+ sinha; -^{k' ■ v - k'^v^) - -rrik • iT - /czi^z) ) ipi{k)(fi{k') 
\ kv k V ' 



(59) 



and 



k± 



AttN^R^ZM r dk±dip 

coshw(l - (cosho; - 1)^)A^ ^ (27r)3 

k k' 
--^(po{k')(fi{k) - -^(pi{k')(po{k) 

+(coshu; - 1) (^{vxk- v - kxv'^)(fo{k')'Pi{k) 



+ 



1 



k'v"^ 



[vxk' ■ V - k'^v'^)(po{k)(pi{k') 



' k ^ k' ^ 

+ sinhu; ( -^{k' ■ v - k'^v^) - -r^{k ■ v - kr^v^) ] fi{k)ipi{k') 
\/cf kv ' 



(60) 



A-kNqR^ZM r dk^dLp 

cosha;(l - (coshw - 1)^)A^ J (2vr)3 

cosh a; - 



ki. 



kk' 
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+ sinh u! 



kv k'v 



(61) 



The valence anti-quark corresponding functions are obtained from these by the fohowing transfor- 
mations 

H\{x,i,t) = -n\{-x,i,t). 

(62) 

The calculation thus far suffers from a traditional problem, namely the so called support problem, 
i.e. the parametrizing functions are non-vanishing outside the physical range [—1,1]. For simplicity 
we use a generalization of the prescription of Ref. ( |2^] ) given by 

which limits the functions to the adequate interval [—1,1], but does not avoid that the quark (anti- 
quark) contribution extends into the negative (positive) x region: our partons are only quark (anti- 
quark) valence partons and therefore the functions should not extend to negative (positive) x. However 
this deffect has a minor impact on the final result. 

Using the above expressions adequately modified by the support prescription and saturating the 
spin flavor degrees of freedom of the pion wave functions we obtain the pion parametrizing functions 
within the boosted scheme which are shown in Figs.||, ^, and |^. In them one can see how the region 
around a; = is the most problematic, but does not affect the asymmetries in an important manner. 

For the sake of completeness and complementarity we have performed the calculation also in the 
unboosted Peirls-Yoccoz scheme. The latter has no support problem, but lacks recoil corrections. 
In this case we obtain the following equations, 

nAx,tt) - NiAnR) \'^\j2^i W) l^2(A/2)P |(l-x)| 

/ %y 
X \ipQ{k)ipQ{k') + —ipi{k)ipi{k') 

- 'y\ oik')Mk) - ^^y\ o{k)Mk')y (64) 
^ 2 ^ ^ Jo ' ^'(27r)2io (27r) |02(A/2)|2 \{1 - x) 



X ll^Q(k')Mk) + l^Mk)Mk')], (65) 



and 



nAx,^,t) - -N{A.R)Mj^ 1^^1(2^70 W) I02(A/2)P |(l-x)| 

X ^^,{k)^,{k'). (66) 



11 



Here the normalization of the wave functions reads: 

, , , ,,n 47ri?3 1 r dv ,2uv 

m{v)\ 



12 • 2 ^2 / sm(— )r2(t;), 

(W^ — SlW^LOy U Jo V ijj 



l(p) 



with 



vdv sin{^) 
(u;^ — sin?u:) Jo uP' 



and the funtion T{v) given by 



T(v) 



UJ 



sin^ UJ 



UJ 



sin 2v 



1 sin 2uj 

2 ^ 2uj 



1 sin 2uj sin uj 



^^^'"+2+ 2UJ 



-V 



UJ^ 



(67) 
(68) 

(69) 
(70) 



The results corresponding to the Peirls-Yokkoz scheme are shown in Figs. ^ |9|, and 
For t — > both approaches are almost the same. However as t grows they become different, and 
in particular the H3 parametrizing function, which is strongly dependent on the boost, becomes very 
small in the boosting scheme. This is the reason behind the smallness of the twist-3 contribution to 
the asymmetry in the latter. 

As a final check we calculate the hadron sum rules that arise from T— invariance [15, [^], i.e. 



-1 



{x,0&{-^ < X < 1) + -7^f")(x,0e(-l < X < e) = F^{t) 



dx{ n'l^i^{x,oei-C < X < 1) - w;:f^(x,oe(-i < x < o 



(71) 



(72) 



We reproduce these sum rules with good precision within our model calculations. In particular when 
t — > we get for the pion form factor ( [n]) 0.94 instead of one. 



6 Concluding Remarks 

The study of DVSC has shown that electromagnetic gauge invariance requires twist-3 contributions. In 
order to check this result experimentally we have analyzed the scattering of linearly polarized electrons 
off hadrons and demonstrated that the SSA is in principle sensitive to the twist-3 contribution. In 
order to be quantitative we have been forced to calculate GPDs, in particular, certain parametrizing 
functions which characterize the needed GPDs. To do so we have performed a calculation of the 
required lightcone matrix elements in the MIT bag model with both boosted and unboosted wave 
functions. 

In Figs.^ — |l^ we show all the parametrizing functions. We have performed the calculation in two 
complementary schemes. The boosting scheme, which takes proper care of the recoil of the pion, but 
does not deal in an exact manner with the center of mass problem and the Peirls Yoccoz scheme, 
which contains no treatment of the recoil, but deals adequately, for small t, with the center of mass 
problem. Certainly the twist-2 Hi is the largest, but the twist-3 ones are non negligible and therefore, 
with an appropriate choice of kinematics, they could be even dominant. However we have restricted 
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our choice to the kinematics which can be achieved in HERMES. For the nucleon, which contains only 
valence quarks in our scheme, the contribution of the twist-3 parametrizing functions will be larger. 

In Figs. and § we show different aspects of our study of the SSA. Fig.^ shows the SSA obtained 
by taking all contributions into account. It is small but measurable with todays high luminosity beams 
and efficient detectors. In Fig.|3| and ^ we depict the pure twist-3 contribution to the SSA, which is 
not small, 15 % at the peak within the unboosted scheme and is much less, 1 % at the peak within 
the boosted scheme. Thus the implementation of gauge invariance is not only theoretically relevant 
but also quantitatively rather important. 

Our model for calculating the parametrizing functions lacks at present one crucial ingredient, 
namely Renormalization Group Evolution. The necessary ingredients to implement such a program to 
next to leading order are not available. Results, however, might be strongly affected by evolution as 
has been the case for other structure functions [28|. Thus until the actual calculation has been evolved 
to the energy regime under scrutiny we will not be fully certain about the experimental relevance of 
this observable. 

The future of the present calculation resides in generalizing our results to the nucleon. For un- 
polarized nucleons, the only difference in the treatment, arises from the kinematics which is more 
complicated because we cannot neglect the nucleon mass. This fact, however, should not change the 
qualitative features of the present results dramatically. For the reasons mentioned above, we expect 
larger values both for the twist-2 contributions, because there are more scatterers, as for the twist-3 
ones, because the contribution of the Ha parametrizing function will not be reduced. The next step 
should be to proceed to the study of polarized nucleons. In this case the number of observables 
increases dramatically due to the spin structure of the target and a study is under way aiming to 



separate them in concrete experiments and estimate their values [29|. 

Our present work shows once again that factorization allows the use of models and perturba- 
tive QCD in a consistent fashion generating a predictive scheme which is useful in guiding future 
experimental developments. 
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Figure 2: The SSA parameter as function of angle (j). The kinematic region of the calculation is defined 
by: t=-0.1 GeV^, x=0.3, ^=0.3, S=22. GeV^. Solid curve: the result within the unboosted scheme. 
Dashed curve: the resilt within the boosted scheme. 
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Figure 3: The unboosted scheme SSA parameter from only twist-3 contribution as a function of (f) for 
the same kinematics as in Fig.0. 




Figure 4: The boosted scheme SSA parameter from only twist-3 contribution as a function of for 
the same kinematics as in Fig.^. 
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Figure 5: The boosted scheme generahzed parton distribution Hi for the vr^ case at t 
and ^ = 0.3. 
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Figure 6: The boosted scheme generahzed parton distribution H'^ for the vr^ case and for the same 
kinematics as in Fig]5|. 
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Figure 7: The boosted scheme generahzed parton distribution Ha for the vr^ case and for the same 
kinematics as in FigjBj. 
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Figure 8: The unboosted scheme generalized parton distribution Hi for the vr^ case at t 
and ^ = 0.3. 
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Figure 9: The unboosted scheme generahzed parton distribution for the vr^ case and for the same 
kinematics as in FigM 
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Figure 10: The unboosted scheme generahzed parton distribution Ha for the vr^ case and for the same 
kinematics as in the Figj8|. 
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